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Abstract 

Monte Carlo (MC) simulations and series expansion (SE) data for the energy, 
specific heat, magnetization and susceptibility of the ferromagnetic 4-state Potts 
model on the square lattice are analyzed in a vicinity of the critical point in order 
to estimate universal combinations of critical amplitudes. The quality of the fits 
is improved using predictions of the renormalization group (RG) approach and of 
conformal invariance, and restricting the data within an appropriate temperature 
window. 

The RG predictions on the cancelation of the logarithmic corrections in the 
universal amplitude ratios are tested. A direct calculation of the effective ratio of the 
energy amplitudes using duality relations explicitly demonstrates this cancelation 
of logarithms, thus supporting the predictions of RG. 

We emphasize the role of corrections and of background terms on the determni- 
nation of the amplitudes. The ratios of the critical amplitudes of the susceptibilities 
obtained in our analysis differ significantly from those predicted theoretically and 
supported by earlier SE and MC analysis. This disagreement might signal that the 
two-kink approximation used in the analytical estimates is not sufficient to describe 
with fair accuracy the amplitudes of the 4-state model. 
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1 Introduction 

In a first paper [1] (hereafter denoted as I), we studied the universal combinations of 
critical amplitudes of the 3-state Potts model. The present paper is devoted to a similar 
analysis in the 4-state case, which is much more involved due to the presence of logarithmic 
corrections strongly influencing the critical behavior. 

We analyze numerical data obtained in Monte Carlo (MC) simulations using the 
Wolff [2] single-cluster algorithm and also the series expansion (SE) data available in 
the literature. In the following, we refer to the data type as MC and SE data, respec- 
tively. For comparison with our own results, we shall also reconsider the data obtained 
in MC simulations using the Swendsen-Wang cluster algorithm [3] by Caselle, Tateo, and 
Vinci [4] and indicated as CTV data. 

Our motivation in using different data sets is to achieve a better control of the critical 
behavior, since one may expect, for the three sources, some differences in the critical region 
due to the different interplay of the finite size effects. In addition, one can apply different 
techniques to the data analysis: the fits in the case of the MC data and the approximant 
technique in the case of the SE data. The consistency of the final results will increase 
our confidence. We care so much because the presence of logarithmic corrections makes 
the numerical determination of the critical behavior of the 4-state Potts model a rather 
delicate task. 

To be even safer, two different approaches are used, which were successfully applied 
also to the q = 3 state Potts model in I. First we estimate the critical amplitudes, which 
are then used to compute universal ratios. Second, besides the direct determination of 
the amplitudes themselves, we estimate ratios of critical amplitudes, constructing effective 
ratio functions, and computing their limiting values at the critical point. This provides 
a direct estimate of universal ratios. Analyzing the renormalization group equations, we 
have shown in an Appendix (see also [5,6]) that, in the absence of any regular background 
term, the logarithmic corrections cancel in the effective ratio functions. 

The Hamiltonian of the ferromagnetic Potts model [7] reads as 



where Sj is a "spin" variable taking integer values between and q — 1, and the sum is 
restricted to the nearest neighbor sites (ij) on a lattice of sites with periodic boundary 
conditions. The partition function Z is defined by 



in zero magnetic field, the model is self-dual. Denoting by (3* the dual of the inverse 
temperature /?, the duality relation 




(1) 




conf 



with P = l/fc^T^, and the Boltzmann constant (fixed to unity). On the square lattice 



(e/^ _ 1) - 1) = g 



(3) 



^According to the usual terminology, the inverse temperature and the critical exponent of the magne- 
tization are denoted by the same symbol /?, since there is no risk of confusion in this context. 
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determines the critical value of the inverse temperature [7] (3c = ln(l + ^yq) ~ 1.09861. 
Dual reduced temperatures r and t* can be defined by 

/3 = /3,(l-r) and (3* = (3^(1 + r*). (4) 
Close to the critical point, r and r* coincide through first order, since r* = t + ^^^^'^^\ ^ + 

The critical amplitudes and the critical exponents describe the singular behavior of 
the thermodynamic quantities close to the critical point. For example the magnetization 
M, the (reduced) susceptibility x and the specific heat C of a spin system in zero external 
fielci behave asn 

M(r) = B{-tY (1 + corr. terms) , r < (5) 
xij)± = r-|-|r|~'^ (1 + corr. terms) , (6) 

C(r)± = — |r|^" (1 + corr. terms) , (7) 



a 

Here r is the reduced temperature t = [T — Tc)/T and the labels ± refer to the high- 
temperature (HT) and low-temperature (LT) sides of the critical temperature Tc. In 
addition to the mentioned observables, for the Potts models with g > 2 a transverse 
susceptibility can be defined in the low-temperature phas^ 

Xt(t) = ^t{—'t)~'^ (1 + corr. terms) . (8) 
The critical exponents are known exactly for the 2D Potts model [8-12]: 

where y is related to the number of states q of the Potts variable by 

cos— = -v/g. (10) 

The standard exponents follow from = {1 — a)/!/ and = 13/v. The central charge of 
the corresponding conformal field theory is also simply expressed [11, 12] in terms of y 

c^,-^. (11) 



Analytical estimates of critical amplitude ratios for the g-state Potts models with 
g = 1, 2, 3, and 4 were recently obtained by Delfino and Cardy [13]. They used the exact 



^In this paper we only deal with the physical properties in zero magnetic field. 

■^Note that for simplicity we have dropped for the moment the multiplicative logarithmic corrections 
and allowed only for additive corrections. 

^In the following we shall use the notations Tl or Vt for the longitudinal or transverse susceptibility 
amplitudes in the low-temperature phase. When still used, r_ is identified with F^. 
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2D scattering field theory of Cliim and Zamolodcliikov [14] and estimated tlie ratios using 
a two-kink approximation for 1 < g < 3. For 3 < g < 4, tfiey considered botli tlie two-kink 
approximation and tlie contribution from the bound state. For g = 4 this approximation 
leads to the value c = 0.985 for the central charge, to be compared to the exact value 
c = 1. Using this approximate value, one can calculate the scaling dimensions from ffTTl) 
and ([9]) obtaining the values Xo- = 0.117 and = 0.577, to be compared with the exact 
values 1/8 and 1/2 respectively. So, the deviation from the exact values becomes as large 
as 6 and 15 per cent, emphasizing the difficulty of the g = 4 case. In the 3-state case the 
situation is much better (see I). 

Let us recall that the existence of logarithmic corrections to scaling in the 4-state Potts 
model was pointed out in the pioneering works of Cardy, Nauenberg and Scalapino [15,16], 
where a set of non-linear RG equations was solved. Their discussion was recently extended 
by Salas and Sokal [17]. 

The universal susceptibility amplitude ratio Tj^/Tl was calculated by Delfino and 
Cardy in [13] for both the 3-state and 4-state Potts models. Later Delfino, et al. [18] 
estimated analytically also the ratio of the transverse to the longitudinal susceptibility 
amplitudes Tt/Tl. The values obtained in the 4-state case are 



In this latter paper [18], the results of MC simulations were also reported, but they 
were considered inconclusive by the authors. More recently Delfino and Grinza report 
compatible values in the case of the Ashkin- Teller model, using the same technique at the 
same level of approximation [19]. 

Another contribution to the study of the amplitude ratios in the 2D 4-state Potts 
model was reported by Caselle, et al [4]. These authors presented a MC determination 
of various amplitudes. In particular, their estimate of the susceptibility amplitude ratio 
r^/Fi = 3.14(70) is in reasonable agreement with the theoretical estimate of Delfino and 
Cardy, in spite of a somewhat controversial [18] use of the logarithmic corrections in the 
fitting procedure. 

Enting and Guttmann [20] also analyzed SE data for the 4-state Potts model and 
found 



results which are compatible with the predictions of [13] and [18]. Their series analysis 
does not rely on differential approximants, but, in the hope to achieve better control of 
the log-corrections of the g = 4 case, they address directly the asymptotic behavior of the 
series coefficients. 

In the present paper we present more accurate MC data supplemented by a reanalysis 
of the extended series derived by Enting and Guttmann [20]. We address the following 
question: Is it possible to estimate the infiuence of the logarithmic corrections on the fit 
procedure? Is it possible to devise some procedure in which the role of the logarithmic 
corrections is properly taken into account? 



F+/Fi = 4.013, Tt/Tl = 0.129. 



(12) 



F+/Fi = 3.5(4), Tt/Tl = 0.11(4), 



(13) 
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In the rest of the paper, we shall be concerned with the following universal combina- 
tions of critical amplitudes 

Tl' Tl' ^" 52 ' ^c- ^2 

where the last two are a consequence of the scaling relatior@ a = 2 — 2/3 — 7. To the 
various critical amplitudes of interest, A±, T±,. . . , we have associated appropriately de- 
fined "effective amplitudes", namely temperature-dependent quantities A±{t), T±{t), . . ., 
which take as limiting values, when r — > 0, the critical amplitudes A±, r±, . . . . In order to 
avoid any risk of confusion between the critical amplitude and the corresponding effective 
temperature-dependent amplitude, reference to these temperature-dependent quantities 
is always made with their explicit r— dependence. 

Also in this paper, we make use of the duality relation in order to improve the estimates 
of the ratios between effective amplitudes measured at dual temperatures. In the case of 
the 4-state Potts model, this procedure would even eliminate all logarithmic corrections 
from the fit in absence of background contributions, which unfortunately do exist for 
most quantities! We again use the duality relation to estimate the correction-to-scaling 
amplitudes in the behavior of the specific heat and of the susceptibility. For this purpose, 
we compute ratios also on the duality line, e.g. the susceptibihty effective-amplitude ratio 
T+{t)/Tl{t*) = x+{P) / Xl{P*) as the ratio of x+iP): the high-temperature susceptibility 
at inverse temperature P, and of Xl{P*), the low-temperature susceptibility at the dual 
inverse temperature (3*. Furthermore, we show analytically that the leading logarithmic 
corrections cancel on the duality line for the ratio of the specific-heat amplitudes as 
extracted from the energy at dual temperatures. 

As a final result of our analysi^, we propose estimates of the susceptibility critical- 
amplitude ratios T^/Tl = 6.49(44) and Tt/Tl = 0.154(12) which are significantly differ- 
ent both from the predictions eq. (fT2l) of [13,18], and from the numerical estimates of [4,20]. 
The deviation from the numerical estimates of other authors might be explained by the 
complicated logarithmic corrections used to fit the data in [4,20]. The difference from the 
theoretical predictions might be due to the limited accuracy of the approximation scheme 
used for g > 3. 

In conclusion, obtaining an accurate (say, within a few per cent) and generally accepted 
approximation of the critical amplitude combinations for the four-state Potts model still 
remains an open issue both theoretically and numerically. 



2 Computational procedures 
2.1 Monte Carlo simulations 

We use the single-cluster Wolff algorithm [2] for studying square lattices of linear size 
L with periodic boundary conditions. Starting from an ordered state, we let the system 

^We refer the reader to the review Ref. [21] for a detailed discussion of the universality of the critical 
amplitudes ratios. 

^The figures given here are an average between the MC and the SE determinations [5]. 
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equilibrate in 10^ steps measured by the number of flipped Wolff clusters. The averages are 
computed over 10^ — 10^ steps. The random numbers are produced by an exclusive-XOR 
combination of two shift-register generators with the taps (9689,471) and (4423,1393), 
which are known [22] to be safe for the Wolff algorithm. 

The order parameter of a microstate M(t) is evaluated during the simulations as 

H=«"^"/^-', (15) 
q-1 ' ^ ' 

where N„i is the number of sites i with = m at the time t of the simulation [23], and 
m G [0, 1, (g — 1)] is the spin value of the majority of the spins. N — L"^ is the total 
number of spins. The thermal average is denoted M = (M). 

Thus, the longitudinal susceptibility in the low-temperature phase is measured by the 
fluctuation of the majority spin orientation 

ksT XL = ^{{N'J - {NJ') (16) 

and the transverse susceptibility is deflned in terms of the fluctuations of the minority of 
the spins 

while in the high-temperature phase x+ is given by the fluctuations in all q states, 

kBTx+ = ^Y.^{Ni)-{N,)'), (18) 

/i=0 

where N^^ is the number of sites with the spin in the state /i. Properly allowing for 
the finite-size effects, this definition of the susceptibilities is, in both phases, completely 
consistent with the available SE data [24]. 

The internal energy density of a microstate is calculated as 



^--^E^^^s, (19) 



N 

{ij) 

its ensemble average is denoted by = (E) and the reduced specific heat per spin measures 
the energy fluctuations, 

BE 

{ksTf C^-— = (E^) - (E)^. (20) 

We have simulated the model on square lattices with linear sizes L = 20, 40, 60, 80, 
100 and 200. In each case, we have measured the physical quantities within a range of 
reduced temperatures called the "critical window" and deflned as follows. Assuming a 
proportionality factor of order 1 in the definition of the correlation length, the relation 
L < ^ oc |r|^^ yields the value of the reduced temperature at which the correlation length 
becomes comparable with the system size L and thus below which the finite-size effects 
are not negligible. This value defines the lower end of the critical window and avoids finite 
size effects which would make our analysis more complex. The upper limit of the critical 
window is fixed for convenience at r = 0.20 — 0.25. 



version of September 26, 2008, Nancy 



7 



2.2 Series expansions 

Our MC study of the critical amplitudes will be supplemented by an analysis of the 
high-temperature and low-temperature expansions for g = 4, recently extended through 
remarkably high orders by Briggs, Enting and Guttmann [20,25]. In terms of these series, 
we can compute the effective critical amplitudes for the susceptibilities, the specific heat 
and the magnetization and extrapolate them by the standard resummation techniques, 
namely simple Pade approximants (PA) and differential approximants (DA), properly 
biased with the exactly known critical temperatures and critical exponents. 

The LT expansions are expressed in terms of the variable z = exp(— /5). In the q = 4 
case, the expansion of the energy extends through z^^. For the longitudinal susceptibility 
the expansion extends through z^^, and for the transverse susceptibility through z*^^. In 
the case of the magnetization, the expansion extends through z^^. The HT expansions 
are computed in terms of the variable v = [1 — z)/{l + {q — l)z). They extend up to f 
in the case of the energy and up to f for the susceptibility. 

It is useful to point out that, for convenience, in Ref. [20] the product of the suscep- 
tibility by the factor g^/(g — 1), rather than the susceptibility itself, is tabulated at HT, 
because it has integer expansion coefficients. For the same reason, at LT the magneti- 
zation times q/{q — 1) is tabulated. Therefore the appropriate normalization should be 
restored in order that the series yield amplitudes consistent with the MC results. 

As a general remark on our series analysis, we may point out that the accuracy of 
the amphtude estimates given in Ref. I for the q = 3 case is good due to the relatively 
harmless nature of the power-like corrections to scaling, while in the g = 4 case the 
mentioned resummation methods cannot reproduce the expected logarithmic corrections 
to scaling and therefore the extrapolations to the critical point are more uncertain. In 
this case we have tested also a somewhat unconventional use of DA's: in computing the 
effective amplitudes, we only retain DA estimates outside some small vicinity of the critical 
point, where they appear to be stable and reliable. Finally, we perform the extrapolations 
by fitting these data to an asymptotic form which includes logarithmic corrections. We 
shall add further comments on the specific analyses in the next sections. 



3 Critical amplitudes of 4-state Potts model 

3.1 Expected temperature-dependence of the observables 

In the case of the 4-state Potts model, we have y = from (ITU]) and the second thermal 
exponent [11,12,26] y^p^ = — 4?//3(l— ?/) vanishes. Accordingly, the leading power-behavior 
of the specific heat (and of other physical quantities) is modified [15] by a logarithmic 
factor 

C{r) = — |r|-"(- In \r\)-' X^^„.{- In |r|) + y,t{\r\). (21) 
a 

The exponent a takes the value 2/3 and A'corr(— In |t|) contains correction terms with 
powers of |r|, — In |r| and ln(— In |r|). It may also contain non-integer power corrections 
due to the higher thermal exponents [8-12] y^j^^ or to other irrelevant fields, as well as 



version of September 26, 2008, Nancy 



8 



power corrections due to the identity field. 3^6* (|t|) contains all regular contributions and 
is referred to as the background term. 

Extending the pioneering work of Cardy, Nauenberg and Scalapino [15, 16], Salas and 
Sokal [17] obtained a set of non-linear RG equations. In the appendix (see also Refs. [5, 
6]), we derive from this set of equations a closed expression for the leading logarithmic 
corrections, which is more suitable to describe the temperature range accessible in a 
numerical study, than the asymptotic form given by Salas and Sokal, 



a 



ln|r|) 



-1 



31n(-ln|r| 
~ 2 -Inlrl 



O 



In Irl 



(22) 



which is the first term of a slowly convergent expansion of Xcorr{ — In |t|) in logs. We have 
observed that the following expansion (see Appendix) is better behaved in the temperature 
window near the critical point accessible by MC and SEs 



S(- 



C{r) 
ln|r|) 

.ln|r|) 



A 



— |r|-2/3e;-i(-ln|r|), 
a 

ln|r|) X ^(-ln|r|) x J^(-ln|r|), 
31n(-ln|r|)' 



1 + 



4 



■ In Irl 



^_3M 

4 -ln|r 



ln|r|) 



ln|r|) 



^(-ln|r|) 



1 



Ci C2ln(-ln|r|) 



-In 



(-ln|r|)^ 



-1 



(23) 
(24) 



(25) 
(26) 



The function £^(— ln|r|) contains the exact form of the leading terms with universal co- 
efficients predicted by RG. The remaining part is made of log terms, whose coefficients 
involve the non-universal dilution field ipQ. The multiplicative function JF(— ln|r|) mim- 
ics, in a given temperature range close to the critical point, the higher-order terms of this 
non- universal part, which is a slowly convergent series in powers of logs starting with an 

O ^i^q^i j term. The values of Ci and C2 in (!26|) should not be considered as "real" ampli- 
tudes of correction-to-scaling terms, but only as "effective" parameters. A "zeroth-order" 
analysis may be performed taking Ci = C2 = 0, i.e. T{— In |r|) = 1. A more refined esti- 
mate follows from the analysis of the magnetization, as explained in the section [3^2] and in 
the Appendix. The absence of a constant background term is a simplifying feature in the 
analysis of the magnetizatioE0. The difference in behavior between the two types of ex- 
pressions ( !22l) and (|23l) for the specific heat is illustrated by the fact that, for example, at 

0.10, 

0.786 



a typical value of temperature numerically accessible without finite-size effects, r 



we have 1 



3 ln(- 
2 



ln|r|) 



0.457, while 1 



3 ln(- 
4 



ln|r|) 



1 + 



3ln(- ln|r|) 



-ln|r| — w.-nui, v»±±iiv. ^ -ln|T| j \^ ^ i -ln|r| 

and the two functions are not simply proportional to each other in the typical range 



^The effective amplitudes are constructed by dividing the corresponding physical quantity by the lead- 
ing terms (main power dependence with known logarithmic corrections), and, therefore the background 
terms, if present, will be divided by the logarithmic terms. These terms may seriously complicate the 
analysis of the limiting behavior of the effective amplitudes. 
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r = 0.02 — 0.20. Fitting the numerical data with one or the other choice may thus spoil 
the outcome for the leading amplitude. 

Similar expressions of the logarithmic corrections are obtained for the other physical 
quantities in Eqs. (17911851) of the Appendix. 

Further corrections to scaling may also be present in Eq. fl23|) . They are discussed 
in the Appendix and may be of the form 02/3 |Tp/^, ah^\T\^'^3 or 0(^,3 |r| ^'''3, as well as 
powers of these terms, where =3/4 and A^,^ = 5/3. Pure power corrections and 
background terms may also be needed. Here we also stress that the inclusion of a leading 
correction in a|rp/^ and of analytic terms seems to be necessary according to the papers 
by Joyce [27,28], where the magnetization of a model, expected to belong to the 4-state 
Potts model universality class, is shown to have an expression of the form 

M(-|r|) = |r|^/^2(/o(r) + 1x^/^1^ + 1x^/^2^) (27) 

with fi{T) analytic functions when r — > 0~. The correction exponents are obtained from 
the table of the conformal scaling dimensions by Dotsenko and Fateev [11,12], but not all 
of them are necessarily present. However, at least the presence of the exponents 2/3, and 
possibly of 4/3, seems to be needed in order to account for the numerical results. Caselle 
et al. [4] also considered an 02/3 1'rp''^ term to fit the magnetization. The parameter 
is a priori possibly needed only for magnetic quantities, while the corrections in a^,_^ will 
systematically be dropped in our fits, since they are sub- sub- dominant. 

In conclusion, the most general expression that we will consider is the following: 

Obs.(±|r|) ^ Ampl. x |r|"^ x ^*(— In |r|) x (1 + corr. terms) + 

+ backgr. terms, (28) 

corr. terms = a2/3|rp/^ + 6-t|r| + a4/3|r|^/^ + . . . , (29) 

backgr. terms = Dq + Di\t\ + . . . (30) 

where ^(— ln|r|) is defined by Eqs. (12411261) . while < and ^ stand for exponents which 
depend on the observable considered. They are all given in the Appendix. 



3.2 The magnetization amplitude 

The amplitude B of the magnetization is defined by the asymptotic behavior (see the 
appendix for details) 



M{-\t\) = B\T\'^'\-\n\T\)-'^ 



31n(-ln|r| 

1 + T- 



4 (-Inlrl 



^ _ 31n(-ln|r| 



4 (-ln|r| 



X 1 + 



4 (- Inlrl 



^(-ln|r|) 



-1/8 



;i + a|r|2/3 + 6|r| + ...)(31) 



We can extract an effective function Teff{ — In jrl) which mimics the real one JF( — In |r|) 
in the convenient temperature range |r| ~ 0.01 — 0.10. This is done by plotting an effective 
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magnetization amplitude 

5e//(-|r|) = M(-|r|)|r|-Vi2(_in|^|)V8 



. , 3 1n(-ln|r|) 
4 (-ln|r|) 



X 1 



31n(-ln|r|) 
4 (-ln|r|) 



-1 



4(-ln|r|) 



(32) 



which is then fitted to the expression 



B 1 



In |r 



+ 



C2ln(-ln|r|) 
(-ln|r|)2 



1/8 



X (l + a|r|2/3 + 6|r| + ...) (33) 



in which we have also included corrections to scaling. As we have already noticed, the co- 
efficients appearing in the function JF are effective parameters adapted to the temperature 
window considered, therefore the values of Ci and C2 have no special meaning. 

In order to analyze the numerical data and to extract the different coefficients, one 
needs very accurate data. As an illustration, in Figure [1] we compare MC and SE data, 
and MC data from Caselle et al. [4]. 




Figure 1: The magnetization M in the critical window region. Our MC data are represented 
by boxes, the MC data from Ref. [4] by stars and the SE data by a soHd hne. 



The behavior of i?e//(— |t|) is shown in Fig. [2l In tabled], we present a selection of 
our fits of MC data to Eq. (1331) . The ffist column of the table indicates the different 
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0,90 -] ^ , ^ , ^ , ^ , ^ , ^ 1 

0,00 0,05 0,10 0,15 0,20 0,25 0,30 



Figure 2: The effective amplitude of the magnetization M. Insert: The magnetization M as 
function of |t|"^''"^^(— In |r|)~-'^/^. (Our MC data are represented by open circles, the SE data by 
boxes and the fit by a solid line). 



choices of the function J-'{ — In |r|). For each line in the tables, several fits have been tried, 
varying the number of points in the interval |r| G [0.005,0.25] (total number of points 
50) and calculating the y^/d.o.f. for each fit. A reasonable balance has to be found 
between the distance of the points from the critical temperature and their number. It 
appeared that limiting the fit window to 20 data points, i.e. to the interval [0.005,0.1], 
gives the best confidence level. This choice is quite satisfactory, since it corresponds to 
a close vicinity of the critical point. The criterion that we adopted in order to select 
a most convincing fit among all possible fits is the stability of the correction-to-scaling 
amplitudes a and h in equation (13T|) when the temperature window is varied, typically in 
the range |r| G [0, 0.06] to |r| G [0, 0.30]. As long as these numbers have not converged to 
given values, we cannot pretend that it is meaningful to include such corrections in the 
analysis. We have to mention that stability of the correction amplitudes is never reached 
if we stop the correction terms in equation ( 1311) at the leading log, or even at the three 
higher log terms (correction function £^(— In |r|)). Then even the leading amplitude B is 
questionable. The effective function J-'e//(— In |r|) is essential to reach convergence of all 
but the two effective coefficients Ci and C2 which still strongly depend on the temperature 
window. 
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For a given type of fit defined in the first column, Tables [T] and [2] show the parameters 
of the fit minimizing ^''^^ fhe minimization is performed by varying the width and 
position of the temperature window. 

Once our choice is made for Ci and C2, the results of the fit for the amplitude B show 
a remarkable stability. The entry called fit # corresponds to a "zeroth-order" fit in 
which the function J-'{—ln |r|) is taken equal to unity. It is presented for comparison, in 
order to emphasize the improvement occurring in the following lines. Fit #1 keeps all 
correction coefficients, while the lines which follow are obtained decreasing the number of 
fit parameters. We favor fits # 2 and #3, where the coefficient b (the amplitude of the 
linear term |r|) is fixed to zero, since from the first fit (#1) we see that leaving b free leads 
to a value close to zero and we consider more reliable a fit with fewer parameters. Fit # 4 
corresponds to another extreme choice including no irrelevant correction at all. It shows 
that the effective function ln|r|) is more important than the irrelevant corrections 
to scaling in order to achieve a stable amplitude B (comparable to the outcomes of fits 
# 1 to # 3). 



Table 1: Various fits of our MC data for the magnetization effective amplitude 
M(-|r|)|r|-^/^2(-_ij^|^|)i/8£-i/8(|_lj^|^|) to the expression Eq. The stars in the first col- 

umn indicate our favorite fits (the reasons for this choice are given in the text). 



Fit # Amplitude Correction terms 



B <x-r-n oc '"'t'"!"!' oc|r|2/3 oc|r| 

— m |t| (— In \t\)^ ' ' ' ' 

MC # 1.1355(4) = - -0.41(1) 0.41(21) 

MC # 1 1.1566(5) -0.740(2) -0.630(51) -0.172(4) -0.018(6) 

MC # 2* 1.1570(1) -0.757(1) -0.522(11) -0.191(2) 

MC # 3* 1.1559(12) -0.88(5) - -0.21(1) 

MC # 4 1.1593(1) -0.25(2) -3.04(5) - - 



Table 2: Same as table [T] for the magnetization data obtained by the SE method. A star in 
the first column indicates our favorite fit with the coefficients Ci and C2 (now fixed and shown 
in bold face in this table and in the forthcoming tables) obtained from the fits of the MC data 
reported in table [H 



Fit # Amplitude Correction terms 





B 


— In T I 


ln(-ln|r|) 
(-ln|r|)^ 


oc |t|2/3 


oc |t| 


SE # 


1.1364(4) 






-0.435(5) 


0.097(10) 


SE # 1 


1.1597(1) 


-0.637(5) 


-1.417(28) 


-0.115(2) 


0.013(2) 


SE # 2 


1.1589(1) 


-0.660(5) 


-1.246(25) 


-0.126(2) 




SE # 2* 


1.1575(1) 


-0.757 


-0.522 


-0.194(1) 




SE # 3 


1.1583(8) 


-0.981(34) 




-0.185(10) 




SE # 3* 


1.1575(1) 


-0.88 




-0.225(1) 




SE # 4 


1.1573(1) 


0.164(7) 


-4.106(18) 
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Table 3: Same as table [T]but using the MC data of Caselle et al. [4] for the magnetization. The 
fit marked by a star is performed using the same coefficient Ci (shown in bold) as in table [TJ 

Fit # Amplitude Correction terms 



B oc -fJ^ (X 1,1) (X |r|2/3 oc Irl 

— ln|T| (_ln|T|)'= II II 

CTV # 1.1386(2) = = -0.518(6) 0.30(2) 

CTV # 1 1.196 ± 1.499 -7.5 ±163 19.6 ± 429 -1.09 ±17.3 1.2 ± 29.2 

CTV # 2 1.148(6) 0.45(56) -3.6 ±1.1 -0.17(4) 

CTV # 3 1.162(9) -1.13(38) - -0.14(12) 

CTV # 3* 1.1561(1) -0.88 - -0.215(1) 

CTV # 4 1.153(2) 0.73(31) -5.5 ± 0.77 - 



The same procedure is now applied to the SE data. It is a non-conventional approach 
to fit SE data, which are usually analyzed by approximant methods, but the exercise 
is tempting. We thus apply exactly the same procedure as in the case of MC data, 
varying the fitting interval and comparing the values of the x^/d.o.f.. The best results 
are collected in table [21 The agreement between the results quoted in the two tables is 
amazing. Not only the amplitudes, but also the correction coefficients are very close to 
each other. In this table (and in the forthcoming tables), we also present the results of 
this analysis of SE data with Ci and C2 fixed to their best values extracted in table [1] 
from the MC data. They are indicated again with a star (SE #2* and SE ^^3*) and in 
order to emphasize the fact that Ci and C2 in this case are not free parameters, they are 
indicated in bold face. A reason for this approach is first to insist on the consistency of 
the results and second to privilege a fit with less free parameters (and in particular no 
free log term). For the magnetization amplitude, a compromise between MC and SE data 
analysis provides our final estimate 

B = 1.157(1). (34) 

As an additional test, we decided to fit also the MC data of Caselle et al. [4], obtained 
using the Swendsen-Wang cluster algorithm, to our functional expression Eq. (l33l) . Notice 
that the various coefficients reported in table [3] (not only the critical amplitude, but also 
the correction terms) are very close to those reported above in table [H In particular this 
fit gives further support to our choice in favor of fits 7^ 2 and 7^ 3. 

In the following tables, we will refer to the best two fits by the labels # 2* and # 
3* (which means that the coefficients Ci and/or C2 are fixed to their values indicated in 
tabled]) when fitting other quantities. Note that the MC data of Caselle et al. [4] are not 
fitted with the choice # 2*, because they do not extend on a range of temperatures wide 
enough to apply the corresponding functional expression. 

3.3 The specific heat and the energy density 

We now turn to the study of the specific-heat amplitudes. Figure [3] shows the effective- 
amplitude ratio A^{t) / A_{—\t\) computed at dual and symmetric temperatures following 
the same prescription as in the case of the 3-state Potts model studied in paper [1]. The 
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same ratio computed from SE data is also shown for comparison. The hnear fit of the 
most accurate MC data (lattice size L = 100 computed with 10"^ measurement steps) 
yields for this ratio the value 0.9999(1), which is remarkably close to the exact value 1 
derived from duality. 




Figure 3: The specific- heat effective-amplitude ratio computed from the energy ratio 
^ J^^^^l^^"!^, as a function of the reduced temperature r on the dual line. The SE data 
are represented by boxes, our MC data (L = 100) by circles. The solid line is given by Eq. (I36p . 
The ratio of the effective amplitudes A^{t) /A-{—\t\) is shown for different sizes in the insert 
(for lattice linear sizes L = 100 (open circles), L = 200 (squares), L = 300 (up triangles), and 
L = 400 (down triangles). The simulations were performed with Nth = 10^ thermalization 
steps and Nmc = 10^ MC steps. The results of a simulation for L = 100 with N^^ = 10^ and 
Nmc = 10^, are represented by closed circles. The SE data are represented by a solid line. The 
dotted line represents a linear fit of the closed circles. 



The ratio 

A+jr) ^ {E{P) - Eo)r"-^ 

A4t*) (Eo-E(/?*))(r*)--i' ^ ' 

where the constant Eq is the value of the energy at the transition temperature, Eq = 
E{j3c) = — 1 — 1/y/q, when expanded close to the transition point leads to 

"^^^^^ - l + la,T + 0{T'+^) (36) 



A_{t*) 3 
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with Oq = —EojSce = '"(-'^v^^ ~ 0.5493. Eq. fl5B]) predicts an asymptotically linear 
r— dependence of the effective-amplitude ratio. This linear dependence is observed in 
Fig. [3] which also confirms that the leading logarithmic corrections in the scaling function 
'^corr(— ln(|T|)) asymptotically cancel in the ratio A^{t)/A_{t*). (Of course this is true 
also for the same quantity computed at symmetric temperatures, since r* ^ |r| for small 
values of r). Figure [3] compares the effective-amplitude ratio obtained from MC simulation 
and SE data with Eq. ( l36i) . The slope of MC data is 1.25, very close to the predicted 
value 1.28. 

Now, we make the natural conjecture (proven in the Appendix in the absence of 
background corrections) that the cancelation of the leading logarithmic corrections will 
also occur for the other ratios. For the leading log- correct ion, — In |r|, and for the next 
correction in ln(— In |r|)/(— In |r|), this can be shown analytically from the RG as first 
indicated by Cardy, Nauenberg and Scalapino in Ref. [15] and by Salas and Sokal in 
Ref. [17]. Our statement is stronger since it extends also to the higher order log-terms, 
such as the next correction in l/(— ln|r|). We believe that equations fl79l) to fl85l) in the 
appendix are exact, and since all the log-terms come from Eq. fl78l) . they should cancel 
in the appropriate ratios (i.e. when the same powers of the dilution field appear in the 
numerator and the denominator. This is always the case when one considers an effective 
combination tending to a universal ratio as r ^ 0). 

According to a RG analysis (see Appendix |A]), we may write the energy in the critical 
region as 





r 


\'/' l + a±| 


r 


2/3 


a{l - a)Pc (- 


-In 


\r\)S{-\n\ 


\rm- 


-In 


\r\) 



In a fixed range of values of the reduced temperature, the "correction function" JF( — In |r|) 
is now fixed and the only remaining freedom is to include background terms (coeff. D) 
and possibly additive corrections to scaling coming from irrelevant scaling fields (coeff. 
a). Therefore, once the function JF(— In |r|) is fixed after our study of the magnetization, 
a reasonable fit of the energy data needs only three parameter^, A±, a±, and D±. The 
next step is the fit of the mean A{t) of the effective amplitudes, 

A{t) = la{l-a)f3,{E4T)~E4-\r\))xg{-\n\r\)/\r\'-'^ 

= ^{E^{r)-E.{-\r\))xg{-\n\r\)/\r\'/\ (38) 

We thus fit the MC data to the expression 

A{t) = A{1 + a|rp/3) + Di\r\ x g{- In |r|)/|r|^/^ (39) 

The log corrections, which indeed cancel in the singular part, unfortunately reappear in 
the background term, albeit suppressed by the power |rp/^. 

We note that A{t) is constructed in Eq. (1551) using the values of energy density com- 
puted at symmetric temperatures. The same quantity constructed from the energy densi- 
ties at dwa/ temperatures (with E_{t*) instead of -E_(— |r|) in Eq. fl38l) ) can also be studied 



^Like in the case of the magnetization, a fourth parameter in fo±|T| appears to be unnecessary. 
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Table 4: Fits of the energy difference (MC data computed at dual temperatures) to the expres- 
sion A{t) ~ Ampl. X (1 + corr. terms) + backgrd. terms x |r|~^/^t/(— In |r|). 



Fit # 


Amplitude 




correetion terms 


background term 




Ampl. 


— In T 


(_ln|T|)^ ^ 1^1 


oc \t\ 


MC # 2* 
MC # 3* 


1.338(3) 
1.316(10) 


-0.757 
-0.88 


-0.522 -4.98(8) 
-4.88(38) 


0.920(13) 
0.899(62) 



and provides better results. In table HI we show our results (the best fit is obtained with 
the choices of Cj's coefficients labeled 2* in tables [1] and [2]). Again, the agreement between 
the two fits is quite good, but this time it is more trivial since the same data set is fitted. 
Taking the average of the parameters from the two fits, we conclude 

A = 1.327(12), (40) 

and a ~ —4.93(23), and Di ~ 0.910(38). Pade approximants of SE data for the specific 
heat provide A = 1.35(1). 

The expansion including corrections to scaling and background terms for the specific 
heat follows from the expressions of the energy density. There is some disagreement 
between our amplitude (^i^a)i3 — 5.922(40) —6.021(13) and the result reported by Caselle 
et al. [4], 6 A ~ 7.80 (36)@. We have to notice that the amplitudes are very sensitive to the 
expression used for the fits. Our choice of effective amplitude in Eq. ( 1551) is supported by 
the quite regular behavior shown in figure HI and also by the natural choice of the fitting 
expression (!39|) . The comparison between the MC data and the resulting fit is shown in 
figure [5l We agree on this point with Enting and Guttmann [20] who emphasized that their 
estimates depend critically on the form assumed for the logarithmic sub-dominant terms, 
and on the further assumption that the other sub-dominant terms, including powers of 
logarithms, powers of logarithms of logarithms etc., can all be neglected. 

3.4 Susceptibilities amplitudes 

3.4.1 High temperature susceptibility amplitude 

We proceed along the same lines as for the other physical quantities and fit the high- 
temperature susceptibility to the expression 

X+ir) = ^+r-''^''g^/\-\nT){l + a+T^/^ + b+T) + D+. (41) 

We can easily obtain the amplitude r+ observing that a single constant as a back- 
ground term is sufficient for the fit (this will also be the case at low temperature). 
The effective amplitude 

re//(r) = X+(T)rV6^-3/^(-lnr) 



Notice that Caselle et al. use a different definition of the energy. 
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Figure 4: The effective amplitude A{t) computed from our MC data using Eq. ([38]) (open 
circles) and from the fit to Eq. ()39p (solid line). 



is represented in Fig. [6] and table [5] collects the coefficients determined by the fits. We 
finally obtain the high-temperature susceptibility amplitude 

r+ = 0.0310(7). (42) 

The value which follows from differential approximants to SE data, although less accurate, 
is consistent with it, r+ = 0.033(2). 

3.4.2 Low temperature susceptibilities amplitudes 

The behavior of the longitudinal susceptibility in the low-temperature phase is less easy 
to analyze [5]. We use the expression 

Xl{-\t\) = TL\Tr'/'g{~ In \T\f\l + ai|r|2/3 + 6^|r|) + D^. (43) 

and the various coefficients are collected in table El For the transverse susceptibility, the 
same procedure leads to the amplitudes also listed in the table El 

One may note that the values of the transverse susceptibility amplitude are more 
stable than those of the longitudinal amplitude, while the estimates of the corrections to 
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Figure 5: The energy differences AE^ and AE^, calculated from our MC data for lattice size 
L = 100 (circles) and the fitted expressions (solid lines). 



scaling are less scattered in the latter case. Our final estimates are 

Tl = 0.00478(24) (44) 

and 

Tt = 0.00074(2). (45) 
DA analysis of SE data gives approximately = 0.005(1). 

3.5 Universal amplitude ratio R~. 

We use the available MC data for C, M, and x fo estimate the universal amplitude ratio 
in the low-temperature phase. To this purpose, we form the function (compare with 
Eq. ([69])) 

p-/ I h 2 Ci-\T\)XLi-\T\) , . 

Rci-M) = ar M^{-\r\) ^ ' 

which is an estimator of the universal amplitude ratio in the limit |r| — 0. As discussed 
in the Appendix, we expect that all sets of logarithmic corrections cancel in this ratio. 
Figure [7] shows with open symbols the combination from Eq. fH6|) for two sets of MC data. 
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those of Caselle et al. (CTV) and our simulations. We may fit these data with correction- 
to-scahng terms starting from |rp/'^ or, assuming in plain analogy with the energy ratio 
that such corrections cancel, with terms starting with |r|. In table [7] we include these fits 
for our MC data set varying the temperature window and the number of correction terms. 

A similar analysis was also performed for the data from Caselle et al (not shown here). 
The results of both analysis are consistent. 

More traditionally, we may evaluate the ratio using the estimated values of the 
amplitudes reported in sections [3.2[ 13.31 and [3.41 (see Eq. (IT^ ). The results will be pre- 
sented later, but anticipating the forthcoming analysis, we can quote as a reliable estimate 
R~ ~ 0.0052 ± 0.0002 (approximants of SE data lead to 0.0050(2)). One can see from 
the table [7] that the ratio obtained from the effective function Eq. (1461) is systematically 
larger. The difference may be explained by the fact that fitting the effective ratio function 
Eq. fH6|) without any logarithmic correction, we assume that the background corrections 
for the longitudinal susceptibility xl and for the specific heat C are small in the criti- 
cal temperature window. While this is indeed the case, these background terms are not 
negligibly small and their presence leads to systematic deviations of the estimates of 
presented in the table [71 The ratio fHBI) may be written as aT'^{Cs + Cbt){Xs + Xbt)/M'^, 
where Cs and Xs are the singular parts of the specific heat and of the susceptibility and 
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Table 5: Fits of the high-temperature susceptibihty to the expression x+i'^) — Ampl. 

^^/^(— Inr) X (1 + corr. terms) + backgr. terms. A star in the first column indicates that the 

coefficients Cj's are those deduced from the MC fits of tabled) 



fit # 


amplitude 




correction terms 




backgr. 




r+ 


oc } 

— Ill T 


In(-lnT) 
(-Inr)^ 


(xr2/3 


ocr" 


MC # 2* 
MC # 3* 
CTV # 3* 
SE # 2* 
SE # 3* 


0.03144(15) 

0.03178(30) 

0.03051(29) 

0.03041(1) 

0.03039(1) 


-0.757 

-0.88 

-0.88 

-0.757 

-0.88 


-0.522 
-0.522 


0.561(60) 

0.53(23) 

1.48(34) 

1.30(1) 

1.67(1) 


-0.053(17) 

0.052(120) 
-0.45(24) 
-0.362(9) 
-0.59(2) 



Table 6: Fits of the low-temperature longitudinal susceptibility to the expression XlI"!''"!) — 
Ampl. X |t|^''/^ X G^^^{— In |r|) x (1 + corr. terms) + backgr. terms and of the low-temperature 
transverse susceptibility to the expression Xt{—\t\) — Ampl. X |t|-7/6 X g3/4(_in|r|) x (1 + 
corr. terms) + backgr. terms. A star in the first column indicates our favorite fit with the 
coefficients deduced from the MC fits of table [TJ 



fit # 


amplitude 




correction terms 




backgr. 






oc } 

— In T 


In(-lnr) 
<^ (-Inr)^ 


ocr2/3 


ocrO 


MC # 2* 
MC # 3* 
CTV # 3* 
SE # 2* 
SE # 3* 


0.00454(2) 
0.00484(3) 
0.00494(3) 
0.00483(1) 
0.00493(1) 


-0.757 

-0.88 

-0.88 

-0.757 

-0.88 


-0.522 
-0.522 


-2.83(3) 

-3.73(14) 

-4.35(15) 

-3.77(3) 

-4.18(5) 


0.050(2) 

0.13(1) 

0.210(19) 

0.116(3) 

0.178(6) 




Ft 


oc } 
— In T 


In(-lnr) 
^ (-Inr)^ 


ocr2/3 


oc rO 


MC # 2* 
MC # 3* 
SE # 2* 
SE # 3* 


0.00076(1) 
0.00073(1) 
0.00073(1) 
0.00073(1) 


-0.757 
-0.88 
-0.757 
-0.88 


-0.522 
-0.522 


-0.805(34) 
-0.25(13) 
-0.577(14) 
-0.369(15) 


-0.0028(2) 
-0.0050(14) 
-0.00373(15) 
-0.00457(16) 



Cfjt and Xbt are the corresponding background nonsingular terms. Eq. (H6l) may be rewrit- 
ten as ar^CsXs(l + Cht/Cs){l + XbtlXs)/^'^- Thus, the background terms Cu and xu 
contribute when divided by the singular terms (or in other words multiplied by factors 
t'^/^Q and r^/^^~^/^, respectively). Clearly, in the critical region, the first factor has the 
dominant contribution. This "large" term may be eliminated completely if we form a 
quantity equivalent to Eq. fj46|) from the energy difference E^{\t\) — Eq instead of the 
specific heat. 



R^, (— |r|) = a{a — 1)/9cT 



(^-(-1 


\t\) - Eo)xlH 


\r\) 


M2(-| 


\r\) 



(47) 



which is shown in the Figure [7] with closed symbols. The extrapolation at r — ^ 0^ is 
obviously different. The results of the fit of MC data to Eq. (147|) are given in table [HI The 
outcome for the universal combination is now fully consistent with the value 0.0052(2) 
and supports our idea that the specific heat background term spoils the behavior of the 
estimator (H^ . 
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Figure 7: The functions Rq{—\t\) (open symbols) and K^*{—t) (closed symbols) which ap- 
proaches the universal ratio Rq as r ^ 0~. Our MC data are represented by circles and the 
CTV data by stars. 



Again, a similar analysis of CTV data (see table [H]) leads to fully consistent results. 



4 Discussion 

Our final goal is the determination of some universal combinations of amplitudes. This 
can be done either directly from the values of the amplitudes listed in the various tables 
of this paper, or also by extrapolating the effective ratios to r = 0. Let us start with an 
estimate obtained by the second method, and let us concentrate on the most controversial 
amphtude ratios, those of the susceptibilities. As we have just shown in the section on R^, 
this method may lead to systematic deviations if the background terms are not handled 
with care. Later on we shall discuss the ratios of the amplitudes listed in the previous 
tables. 

The ratio Tj^/Tl can be estimated from the ratio of the SE of the susceptibility x+ 
at high-temperature and of the longitudinal susceptibility xl in the low-temperature 
phase. We have again two options to form this ratio, either from quantities computed 
at temperatures symmetric with respect to the critical temperature Tc ± r or at inverse 
temperatures /? and 13* related by the duahty relation Figure [H] shows both the ratios 
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Table 7: Estimates of the critical amplitude ratio Rq from our MC data using different 
fits and varying the temperature window. 



r— window 


amplitude 




correction terms 








(X |t|2/3 


(X |t| 


(X |r|4/3 


(X |t|5/3 


0.01 - 0.29 


0.00651(3) 
0.00685(3) 


-3.04(12) 
-4.60(3) 


-1.24(39) 
3.84(5) 


4.18(32) 






0.00619(1) 
0.00591(2) 




-16.13(11) 
-11.01(10) 


29.25(37) 
12.13(16) 


-14.34(31) 


0.01 - 0.10 


0.00628(7) 
0.00654(3) 
0.00627(3) 


-1.54(58) 

-3.66(9) 

-2.70(2) 


-6.78 ± 2.36 
-1.99(18) 


9.41 ±2.53 






0.00618(3) 
0.00611(1) 




-15.90 ± 1.26 
-13.05(17) 


28.64 ± 5.51 
16.10(35) 


-14.09 ±6.18 


0.01 - 0.046 


0.00642(6) 
0.00637(1) 


-3.09(33) 
-2.83(2) 


0.62(80) 








0.00616(4) 
0.00588(4) 




-14.20(94) 
-6.92(18) 


18.99 ± 2.44 





T^{t)/Tl{—\t\) and T^{t)/Tl{t*), while the straight lines are drawn as a guide for the 
eye. It would be naive to take the linear fit too seriously, otherwise one should conclude 
that the background (non-singular) contribution to the ratio is negligible. The value of the 
universal amplitude ratio T^/Tl obtained from the SE data is approximately T^/Tl — 
6.16(1) and 6.30(1), when using respectively the fits MC #2* and MC #3*. We have 
also analyzed the effective-amplitude ratio from MC data obtained by dividing the high- 
temperature reduced susceptibility by the longitudinal reduced susceptibility computed at 
temperatures related by the duality relation. Neglecting the constant background terms 
in the susceptibilities eliminates all logarithms and makes the fit quite simple, leading to 
a ratio in the range T^/Tl ~ 6.30 — 6.60. On the other hand, if we keep in the fit the 
background terms, the logs reappear and we are lead to Tj^/Tl ^ 6.0 — 6.1. 

Thus we get values which are quite different from the analytical prediction T^/Tl = 
4.013 of Delfino and Cardy [13], as well as from the value 3.5(4) estimated by Enting and 
Guttmann from an analysis of the SE data for the susceptibility in both phases, and from 
the value 3.14(70) estimated by Caselle et al. [4]. 

Let us now estimate the effective- amplitude ratio rT'(— |T|)/rL(— |r|). This ratio, 
shown in figure [HI has been computed both by MC simulation, for various lattice sizes, 
and from SE data. Due to the non-singular correction terms, its behavior is far from 
being linear in r as was the ratio x+{\'^\)/Xl{—\t\) (compare with figure [8]). A possible 
explanation is that there might be some symmetry in the correction-to-scaling amplitudes 
occurring in the asymptotic expansion of x+(|r|) and Xl(— |t|), but not of Xt(— |t|), which 
introduces here a stronger background term still containing the logs. Following again the 
same procedure, we arrive at estimates close to 0.146 when neglecting all logarithmic 
corrections, while we get 0.152 — 0.153 when allowing for these corrections and from the 
analysis of SE data. Finally we obtain Tt/^l = 0.151(3) and 0.148(3) from the fits MC 
#2* and 3*. All these values differ from the analytical prediction Tt/Tl = 0.129 of Delfino 
et al. [18] and from the value 0.11(4) of Ref. [20]. 
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Table 8: Estimates of the critical amplitude ratio R^, (Expr. (H7I) ) from our MC data 
using different fits and varying the temperature window. 



r— window 


amplitude 




correction terms 








(X |r|2/3 


(X |r| 


<x |r|4/3 


(X |t|5/3 


u.ui — ■ u.zy 


u.uuooi i^y ) 
0.00558(3) 


— Z.Ool^O J 

-2.97(3) 


0.57(15) 
1.95(2) 


1.11(12) 






0.00527(8) 
0.00508(1) 




-11.44(16) 
-7.29(9) 


20.98(51) 
7.38(14) 


-11.25(42) 


0.01 - 0.10 


0.00558(4) 
0.00553(1) 
0.00535(2) 


-3.28(36) 

-2.77(5) 

-1.98(2) 


3.66 ± 1.45 

1.58(9) 


2.23 ±1.55 






0.00538(2) 
0.00525(2) 




-16.11(99) 
-9.55(24) 


40.39 ±4.30 
11.82(50) 


31.97 ± 4.81 


0.01 - 0.046 


0.00558(3) 
0.00543(2) 


-3.09(20) 
-2.11(3) 


2.34(48) 








0.00535(3) 
0.00511(4) 




-12.40(74) 
-5.04(19) 


18.95 ± 1.99 





Table 9: Estimates of the critical amplitude ratio i?^, from Caselle et al MC data using 
different fits and varying the temperature window. 



r— window amplitude 




correction terms 




R^ 


(X |r|2/3 


(X |r| 


(X Irl-i/a 


0.0058 - 0.029 0.00548(2) 
0.00543(1) 


-2.59(17) 
-2.14(2) 


1.27(48) 




0.00535(1) 
0.00521(2) 




-13.77(66) 
-5.93(23) 


21.79 ± 1.98 



Let us now extract what we believe are more reliable estimates for the universal com- 
binations of amplitudes by a direct evaluation of the ratios of the numbers presented in 
this paper and collected again in table [THl The universal combinations are presented in 
table [m, together with the corresponding results available in the literature. Averaging 
our different results, we quote the following final estimates: 

T+/Tl = 6.49 ±0.44, (48) 

Tt/Tl = 0.154 ±0.012, (49) 

R-^ = 0.0338 ±0.0009, (50) 

R- = 0.0052 ±0.0002. (51) 

These results clearly confirm the above-mentioned limits of effective ratios. 

The main outcome of this work are the surprisingly high values of the ratios T^/Tl, 
^t/^l and R^ (and the low value for Rq), significantly deviating from the predictions 
of Delfino and Cardy. We emphasize that our results are also supported by a direct 
extrapolation of effective-amplitude ratios for which most corrections to scaling disappear. 
We believe that our fitting procedure is reliable, and since the disagreement with the 
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Figure 8: The ratio of the effective amphtudes r_|_(r)/ri(— |r|) obtained from SE data at 
symmetric temperatures (boxes) and dual temperatures (circles) together with the corresponding 
linear fits (solid lines). 



theoretical calculations can hardly be resolved, we suspect that it might be attributed to 
the approximations made in Ref. [13] in order to predict the susceptibility ratios. Even 
more puzzling is the fact that Delfino and Cardy argue in favor of a higher robustness of 
their results for Fj^/r. than for F_|_/F_, while the disagreement is indisputable in both 
cases. Indeed, in the conclusion of their paper, and in a footnote, Delfino et al [18] 
(p. 533) explain that their results are sensitive to the relative normalization of the order- 
and disorder-operator form-factors, which could be the origin of some troubles for q = 3 
and 4 for the ratios F^/F^ and Rc only. 

As a final argument in favor of our results, we may mention a work of W. Janke and 
one of us (LNS) on the amplitude ratios in the Baxter- Wu model (expected to belong 
to the 4-state Potts model universality class), leading to the estimate F+/F_ ~ 6.9 and 
Rq ~ 0.005 [30]. This result, obtained from an analysis of MC data shows a similar 
discrepancy with Delfino and Cardy's result and suggests that further analysis might still 
be necessary. 
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Figure 9: 4-state Potts model. The ratio of the effective amphtudes rr(— |r|)/ri(— |t|) for the 
4-statc Potts model on square lattices of linear sizes L = 20 (boxes), L = 40 (up triangles), 
L = 60 (down triangles), L = 80 (diamonds), L = 100 (stars) computed with Nmc = 10''' MC 
steps, and L = 100 (closed circles) computed with Nmc = 10^. The solid line represents the SE 
data. 
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Table 10: Critical amplitudes and correction coefficients for the 4-state Potts model. They are 
written in the following format Obs.(ib|T|) ~ Ampl. x |r|"^ x Q^{— In |r|) x (1 + corr. terms) + 
backg. terms. The results of the MC analysis of Ref. [4] are compiled together with our results 
obtained by combining the MC and series expansion (SE) data analysis. 



observable 


amplitude 


correction terms, 


xe*(-ln|r|) 


background terms 


source 






^ r 1 


LA. 1 V 1 


1 1 n 
cx \t\^' 


(X |t| 






L.dS^(d)/a(l - a)fjc. 


— 4.98(8) 




Eo 


3.77(6) 


this paper MC# 2 




1.3lD(9j/a(l — a)pc 


—4.88(38) 




Eq 


3.68(30) 


this paper MCv# .V 


X+{r) 


i + = 0.0223(14) 






0.05(14) 


- 






i + = 0.0dl(5j 






- 


- 


[20] 




i + = 0.0dl44(15) 


0.5o1(d0) 




-0.053(17) 


- 


this paper MC# 2 




r+ = 0.03178(30) 


0.53(23) 




0.052(120) 


- 


this paper MC# 3 




i _l_ = 0.03041(1) 


1.30(1) 




-0.362(9) 




this paper bli^# 2 




■p n mnQO/' i ^ 

i + — U.UoUoyi^l ) 


l.b7(lj 




-0.59(2) 




tills paper otiijf o 


Xi(-|r|) 


Tl = 0.00711(10) 






0.02(1) 




[4] 




= 0.0088(4) 










[20] 




Ti = 0.00454(2) 


-2.83(3) 




0.050(2) 




this paper MC# 2* 




= 0.00484(3) 


-3.73(14) 




0.13(1) 




this paper MC# 3* 




= 0.00483(1) 


-3.77(3) 




0.116(3) 




this paper SE# 2* 




= 0.00493(1) 


-4.18(5) 




0.178(6) 




this paper SE# 3* 


Xt{t) 


Ft = 0.0010(3) 










[20] 




Ft = 0.00076(1) 


-0.805(34) 




-0.0028(2) 




this paper MC# 2* 




Ft = 0.00073(1) 


-0.25(13) 




-0.0050(14) 




this paper MC# 3* 




Tt = 0.00073(1) 


-0.577(14) 




-0.00373(15) 




this paper SE# 2* 




Ft = 0.00073(1) 


-0.369(15) 




-0.00457(16) 




this paper SE# 3* 


M{-\r\) 


B = 1.1621(11) 






0.05(14) 




[4] 




B = 1.1570(1) 


-0.191(2) 








this paper MC# 2* 




B = 1.1559(12) 


-0.210(10) 








this paper MC# 3* 




B = 1.1575(1) 


-0.194(1) 








this paper SE# 2* 




B = 1.1575(1) 


-0.225(1) 








this paper SE# 3* 



A Solution of non-linear RG equations and cancela- 
tion of logarithmic corrections in effective- amplitude 
ratios 

For the 4-state Potts model, the non-linear RG equation for the relevant thermal and 
magnetic fields (p and h, with the corresponding RG eigenvalues and i/h, and the 
marginal dilution field ip are given by 

{y4, + y4>i^'^)(l>, (52) 
{yh + yh'4>i')h, (53) 

g{^) (54) 

where h is the length rescaling factor and / = ln6. The function g{ip) may be Taylor 
expanded, g{il)) = |/^2^^(1 + ^^ip + ...). Accounting for marginality of the dilution field, 

there is no linear term. The first term has been considered by Nauenberg and Scalapino, 
and later by Cardy, Nauenberg and Scalapino. The second term was introduced by 



d0 
dh 

d\n.h 
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Table 11: Universal combinations of the critical amplitudes in the 4-state Potts model. 



A+/A- r+/TL rr/TL Rq Rq source 

T OT3 0l29 0.0204 0.00508 [13, 18] 

- 4.02 0.129 - - [19] 

3.14(70) - 0.021(5) 0.0068(9) [4] 

3.5(4) 0.11(4) - - [20] 



1.000(5) 6.93(6) 0.1674(30) 0.03452(25) 0.00499(3) this paper MC# 2* 

1.000(13) 6.57(10) 0.1508(26) 0.03439(63) 0.00524(9) this paper MC# 3* 

6.30(1) 0.1511(24) 0.03336(9) 0.00530(3) this paper SE# 2* 

- 6.16(1) 0.1481(23) 0.03279(24) 0.00532(5) this paper SE# 3* 



Salas and Sokal. In this appendix, we slightly change the notations of Salas and Sokal, 
keeping however the notation i/ij for all coupling coefficients between the scaling fields i 
and j. These parameters take the values ?/</,^ = 3/(47r), yh^ = l/(167r), y^2 = 1/n and 
Vips = — l/(27r^), while the relevant scaling dimensions are = = 3/2 and yt = 15/8. 

The fixed point is at = = 0. Starting from initial conditions 00, h^, the relevant 
fields grow exponentially with /. The field is analytically related to the temperature, 
so the temperature behavior follows from the renormalization fiow from 0o ~ 1^1 up to 
some = 0(1) outside the critical region. Notice also that the marginal field ip remains 
of order O('0o) and ipo is negative, ipo = In zero magnetic field, under a change of 

length scale, the singular part of the free energy density transforms as 

/(^o,0o) = e"^'/(^,0), (55) 

where D = 2 is the space dimension. Solving Eq. fl52|) leads to ln(0/0o) = y(i>l + y<j>ii) J if^dl 
where the last integral is obtained from Eq. flM|) rewritten as jQipdl = ^ln('0'/'0o) + 
^ lnG(?/'05 ''Z')- Note that G{tlJo,ip) takes the value 1 in Cardy, Nauenberg and Scalapino 
and the value ^v-^+^v-^^" Salas and Sokal. Since this term appears always in the same 

combination, we write z = ipQ/ip, z = ^^^^ and in the same way we set x = 0o/0- One 
thus obtains 

/ = -— Inx + ^^ln^. (56) 

y<t, y^yi,^ 

At the critical temperature = 0, the /—dependence on the magnetic field obeys a similar 
expression and one is led to the equality 

/ = — z/ In X + /i In ^ = —Uh \ny + In z, (57) 

where y = ho/h and for brevity we will denote u = l/y^, = I, /i = = 1 z/^ = l/y^ = 

^ and LLh = ^'"^ = We can thus deduce the following behavior for the free energy 

15 yhy^2 30 ^ ^-^ 

density in zero magnetic field in terms of the thermal and dilution fields, or at the critical 
temperature in terms of magnetic and dilution fields 

/(0o,V^o) = x^'^^-^^/(0,^) (58) 
/(/^o,^o) = y'"''^z-''^'^f{h,^P) (59) 
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The other thermodynamic properties follow by derivation with respect to the scaling 
fields, e.g. E{(f)o,^o) = g|^/(V^o, 0o), or 

i?(0o,V^o) = e-^'^E{<P,^), (60) 
expo 

C(</>o,^o) = W,^), (61) 

dh 

M(0o,^o) = e-^'— M(0,V^), (62) 

X(0o,^o) = (63) 

The derivatives ^ and ^ will thus determine the scaling behavior of all the thermody- 
namic quantities. The first one is obvious, ^ = and for the second one, = y~^, 
we express y in terms of x and z using eq. ( 1571) 0. Altogether, introducing the notations 

^ = yh/y,t> = Z and K = yhy<f>i;/ycf>y4,2 - yh^/y-^,^ = | one has 

/(0o,^o) = x^'^^-^^/(0,V^) (64) 

E(0o,^o) = x^'^-i^-^^E(0,^) (65) 

C(0o,^o) = x^'^-^^-^^C(0,^) (66) 

M(0o,V^o) = x^'^-^ ^-^^+'^M(0, ^) (67) 

X(0o,^o) = x^'^-^^^-^^+2'^x(0,^) (68) 

What appears extremely useful in these expressions is that when defining appropriate 
effective ratio£3, the dependence on the quantity z cancels, due to the scaling relations 
among the critical exponents. This quantity z is precisely the only one where the log 
terms are hidden, and thus we may infer that not only the leading log terms, but all 
the log terms hidden in the dependence on the marginal dilution field disappear in the 
conveniently defined effective ratios. For example in an effective ratio like 

= (69, 

M^{x, z) 

all corrections to scaling coming from the variable z disappear. 

Now we proceed by iterations of / = —ulnx + uln^ and z = ^ i^^t-^^'^^f^'^^?]^ . The 
asymptotic solution of Eq. (15^ i^ 

' ^ + 7^!^H-0(l//)V (70) 



^po 1 - ^oy^2l \ {y^2Y I 



"'""it follows that y — ^^''/^^z'""''/^''-^"^'"''''"''/^'*'^'''^ . 

^^i.e. effective ratios which eventually tend to universal limits when r — > 

^^When y^'i = the asymptotic solution of Eq. ((M)) is simply i_^°y ^i - Thus one can try the ansatz 
■0 — i_^^y (1+^(0) with a small correction X(/) to solve asymptotically the full Eq. ([5^ . Keeping only 
terms of order 0{l~^) at most, we are led to the following expression X'{1) = — {1 / l)X (l) + y^a / {{y^2)'^P) 
where we use X{1) — Y{l)/l to eventually obtain Y{1) ~ {y^^ / {y,p^)'^) InZ. 
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We get for the variable z 

-| _|_ y-tp^ In / 



-ijoy^2u{-\n \t\) 



^ /. In(-lnM) Q 
1/ — ln|r| 



■ln|r| 



1 I ln(-lii|r|) I ^ 

{y^2)^^ -ln|r| V-l°l^l 



(71) 



Similarly, one has the combination 

1 + (y^a/Z/i/^OV^ 1 



1 - 



o 



and eventually one gets for the full correction-to-scaling variable the heavy expression 



(72) 



const X 



const X 



1 + 

ln|r|) 



mM- ln|r|) 



■ln|r| 



y^3 ln(-ln|T|) 



X 



ln|r|)- 



1 + 
1 + 



{y^2rf -in|T| 

3ln(- ln|T|) 



■ln|r| 



3 ln(-ln|T|) 
4 -ln|T| 



(Z/^O^'^I-Inlrl) 

1 



4(-ln|r|) 



X F(-ln|r|) 
x^(-ln|r|) (73) 



where In |r|) is a function of the variable ( — In |r|) only where also appears the non- 
universal constant ipQ. Using Eq. fl671) . we deduce the behavior of the magnetization for 
example 



M(-|r|) = 5|r|i/i2(_in|r|)-i/ 



^ 31n(-ln|r|) 
4 -Inlrl 



1 + 



31n(-ln|r|) 



4 -Inlrl 



1 + 5 



1 



4-ln r 



^(-ln|r|) 



-1/8 



• (74) 



Since all these log expressions are "lazy functions", it is unsafe to expand such terms, 
e.g. ^1 — l^^rii^jTpj ~ 1 + | '"|_~n'||-|^\ since the correction term is not small enough 
in the accessible temperature range |r| ~ 0.05 — 0.10. We can thus only extract an 
effective function jFgjj(— In |r|) which mimics the real one JF(— ln|r|) in the convenient 
temperature range. This is done through a plot of an effective- magnetization amplitude 



i?e//(-|r|) = M(-|r| 



r 



-1/12,' 



^ 31n(-ln|r| 
1 + - ^ , I'l' 
4 -In r 



-1 



which is found to behave as 



ln|r|)i/« 

31n(-ln|r|) 
4 -ln|r| 



Ci C2ln(-ln|r|) 

Tn |r| (— ln|r|)^ 



1 



4-ln 



1/8 



(75) 



(76) 



version of September 26, 2008, Nancy 



30 



from which one deduces that the function JF(— In |r|) takes the approximate expression 



J^(-ln|r|) ~ 1 + 



Ci , C 2ln(-ln|r|) 
;-ln|r|)2 



— In Irl 



+ 



-1 



(77) 



What is remarkable is the stabihty of the fit to Eq. ( 1761) . We obtain (see table [T]) Ci = 
—0.757 and C2 = —0.522 which yields an amplitude B = 1.1570(1). It is also possible to 
try a simpler choice, fixing Ci = and approximating the whole series by the C2— term 
only. We then find the value C2 = —0.88 and this leads to a very close magnetization- 
amplitude B = 1.1559(2). 

For the following, we group all the terms coming from the variable z into a single 
function Q{— In |r|) = (— In |r|) x S{— In |r|) x JF(— In |r|) where 



^(-In 



1 + 



31n(-ln|r|) 
4 -In Irl 



^ 31n(-ln|r|) 
4 -ln|r| 



In r 



(78) 



in terms of which the singular parts of the physical quantities take a very compact form. 



fir 
M{-\t\ 

x+(±kl 

XL{-\r\ 
XT{-\r\ 

E±{±\r\ 
C+(±|r| 



B\ 



1/12^-1/8^ 
.1-7/6 



ln|r|) 
-ln|r|) 



TL\rr^'g'/\-\n\T\) 
rT|r|-^/V/'(-ln|r|) 

a(a — r 



|r|V3^-i(-ln 



a 



ln|r|). 



(79) 
(80) 
(81) 
(82) 
(83) 

(84) 
(85) 



The function S is known exactly while the function needs to be fitted to the numerical 
data. In the same range of values of the reduced temperature, the "correction function" 
^(— ln|r|) is now fixed and the only remaining freedom is to include background terms 
and possibly additive corrections to scaling coming from irrelevant scaling fields Efl. Among 



^•^To introduce corrections to scaling, let us consider the case of an irrelevant scaling field, let say g, 
coupled to the temperature field through 

d4> I , , , T dg 

~J[ =y't>'P + VUVV + y<t,g(P9 and — = ygg 

(A > above plays the role of ~yg/y^, and is thus linked to the corresponding negative RG eigenvalue 
yg). Solving for g gives g = goe^"^ (the irrelevant scaling field decays exponentially when one approaches 
the fixed point). Solving for ij: gives V' 



and for 



I = - ln(0/0o) - ln(^/^o) + ^<?o(e^'' - 1). 

y<p y<ijy-4i^ y4>yg 
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the additive correction terms, we may have those coming from the thermal sector A^^ = 
—py^^, where the RG eigenvalues are y^^ = D — \n? . The first dimension y^^ = y^ is 
the temperature RG eigenvalue, the next one, y^^, vanishes and is responsible for the 
appearance of the logarithmic corrections, so the first irrelevant correction to scaling in 
the thermal sector comes from A^^ = — z^y^^ = 5/3. One can also imagine a coupling 
of the magnetic sector to irrelevant scaling fields. The magnetic scaling dimensions Xa„ 
lead to RG eigenvalues yh„ = D — x„^. The first dimension y^-^ = y^ is the magnetic 
field RG eigenvalue. The second one is still relevant, yh^ = 7/8, and it could lead, if 
admissible by symmetry, to corrections generically governed by the difference of relevant 
eigenvalues {yh^ —yh2)/y<i> = 2/3. The next contribution comes from y^^^ = — 9/8 and leads 
to a Wegner correction-to-scaling exponent [29] A^^ = —pyui =3/4. Eventually, spatial 
inhomogeneities of primary fields (higher order derivatives) bring the extra possibility of 
integer correction exponents yn = —n in the conformal tower of the identity. The first 
one of these irrelevant terms corresponds to a Wegner exponent Ai = — z/(— 1) = 2/3 and 
it is always present. We may thus possibly include the following corrections: 

, the first and third ones being always present, while the other corrections 
depend on the symmetry properties of the observables. 
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